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Abstract. We consider "lattice glass models" in which each site can be occupied by at most one particle, 
and any particle may have at most I occupied nearest neighbors. Using the cavity method for locally tree- 
like lattices, we derive the phase diagram, with a particular focus on the vitreous phase and the highest 
packing limit. We also study the energy landscape via the configurational entropy, and discuss different 
equilibrium glassy phases. Finally, we show that a kinetic freezing, depending on the particular dynamical 
rules chosen for the model, can prevent the equilibrium glass transitions. 

PACS. 64.70.Pf Glass transitions - 64.60.Cn Order-disorder and statistical mechanics of model systems 
- 75.10.Nr Spin-glass and other random models 



1 Introduction 

The thermodynamics of vitreous materials is a long stand- 
ing yet very alive subject of study [1]. In spite of much 
research, it remains unclear whether this "amorphous" 
state of matter can exist in equilibrium; even if it can- 
not, probably the underlying crystalline equilibrium phase 
is irrelevant for understanding the physics of glassy sys- 
tems. As experience has shown in other contexts, a lot of 
understanding can be gained by looking at models sim- 
ple enough to be analyzed but that retain the essential 
physics. Following this strategy, several recent works have 
focused on lattice models for structural glasses. 

Lattice gas models with hard core exclusion, i.e., with 
each site being occupied by at most one particle, designed 
to reproduce the glass phenomenology, fall into two dis- 
tinct classes. A first class consists of kinetically constrained 
models [2] which have a glassy behavior forced by a dy- 
namical constraint on allowed moves, but otherwise trivial 
equilibrium properties. An example is the Kob-Andersen 
model [3] where a particle is allowed to move only if be- 
fore and after its move it has no more than some num- 
' ber m neighboring particles. In this case the slow dy- 
namic displays several remarkable properties like for ex- 
ample an avoided transition toward a cooperative "super- 
Arrhenius" dynamics [4]. Physically, the kinetically con- 
strained models are based on the assumption that the 
glassy behavior is due to an increasing dynamical cor- 
relation length whereas static correlations play no role. 
Thus, the possibility of a thermodynamic glass transition 
is excluded from the beginning. 

In contrast, the models we will discuss here belong 
to a second class where the glassy features are generated 
through a geometric constraint on allowed configurations. 
In this case a thermodynamic equilibrium glass transition 
independent of the chosen local dynamical rules may exist. 



Indeed, as we shall show in the following, it takes place for 
models on a Bethe lattice. Such models were first intro- 
duced in Ref. [5] and some variants have been elaborated 
since then [6,7]. 

In this paper, we study the lattice glass models on 
"Bethe lattices" which are random graphs with a fixed 
connectivity. This kind of approach provides an approxi- 
mation scheme for the lattice glasses on Euclidean lattices 
having the same value of the connectivity, in the same 
way as the Bethe approximation allows one to compute 
approximate phase diagrams of non frustrated systems. 
Furthermore this approximation can be improved system- 
atically, at least in principle, by implementing higher order 
cluster variational methods. But this random graph study 
is also interesting for its own sake. In particular we ana- 
lyze in detail the limit of diverging chemical potential. In 
this limit, one recovers an optimization problem which is 
the lattice version of close-packing of spheres, a problem 
that has challenged mathematicians for many decades [8] 
and is still matter of debate today [9] . 

With respect to their Euclidean counterparts, the lat- 
tice glasses on Bethe lattices have one important differ- 
ence: the existence of the random lattice, even with fixed 
connectivity, forbids crystalline ordering. This loss is also 
an advantage in that the absence of crystal phases makes it 
easier to study the glass phase. Indeed, when the density of 
the system is increased, we find a thermodynamical phase 
transition from a liquid phase to a glass phase. We can 
determine the phase diagram analytically, focusing in par- 
ticular on the liquid to glass transition. In many systems 
this transition is of the "random first order" [10-13] type, 
also called "one-step replica symmetry breaking" [14]. On 
the Bethe lattice there are actually two transitions: when 
increasing the density or chemical potential, one first finds 
a dynamical transition in which ergodicity is broken, then 
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a static phase transition. The intermediate phase is such 
that the thermodynamic properties are those of the liq- 
uid phase, in spite of the ergodicity breakdown. At the 
static transition the entropy and energy are continuous, 
with a jump in specific heat. Therefore the Bcthe-lattice 
glass models provide clear solvable examples of a system 
of interacting particles exhibiting the scenario for the glass 
transition which was proposed in [10-13] by analogy with 
some spin glass systems. This scenario is known to be a 
mean-field one, which does not take into account the mi- 
cleation processes that can occur in Euclidean space. It is 
generally believed that nucleation processes transform the 
dynamical transition into some cross-over of the dynam- 
ics from a fast one to a slow, activated relaxation [10-13, 
15]. Whether the static transition survives in realistic sys- 
tem is unknown so far. In this paper we will not discuss 
the relevance and modifications of the mean field scenario 
when applied to finite dimensional problems since we have 
nothing to add to existing speculations. Let us just men- 
tion that the lattice glasses provide the best examples on 
which these issues can be addressed. The first step of such 
a study is to have a detailed understanding of the mean 
field theory, and this is what we do in the present paper. 

Note that some of the results have appeared in Ref. [5] ; 
here we give a new and extensive discussion on the nature 
of the different possible equilibrium glassy phases. 

On top of the thermodynamic study, we have also stud- 
ied the dynamical arrest of the lattice glasses on Bethe 
lattices. The dynamical arrest depends on the specific dy- 
namical rule that is implemented. We show that this dy- 
namical arrest is in general unrelated to the energy land- 
scape transitions found in the thermodynamic approach. 
In particular, in some models, a dynamical arrest takes 
place at a density smaller than the one of the dynami- 
cal glass transition. In kinetically constrained problems, 
such mean field arrest transitions are known to become 
crossovers in finite dimensional systems. The correspond- 
ing arrest behavior on our lattice glass models is not known 
yet. 

The outline of this paper is as follows. In Sect. 2, we 
introduce the models on lattices of arbitrary types and de- 
fine the relevant thermodynamics quantities. In Sect. 3, we 
address the case of loop-less regular lattices, called Cay- 
ley trees when surface sites are included and Bethe lat- 
tices when surface effects can be neglected. We show that 
low densities correspond to a liquid phase described by 
the Bethe-Peierls approximation, but that inhomogeneous 
phases must be present at higher densities. However, the 
strong boundary dependence of models on Cayley trees 
does not allow one to define such dense phase in the in- 
terior. To overcome this problem, we consider instead in 
Sect. 4 random regular graphs; they share with Cayley 
trees a local tree-like structure but are free of surface ef- 
fects and thus provide a natural generalization of Bethe 
lattices adequate for dense phases. Our study on random 
graphs is performed by means of the cavity method [16] 
which predicts for high densities a glassy phase. In Sect. 5, 
we focus on the close-packing limit and discuss in detail 
the nature of this glassy phase. Finally in Sect. 6, we 



comment on the differences between the equilibrium glass 
transition of our models and the kinetic transitions or dy- 
namical arrests related to specific local dynamical rules. 
In particular we show that for some models the dynamical 
arrest takes place before the equilibrium glass transition. 



2 Lattice Models 

2.1 Constraints on local arrangements 

When packing spheres in three dimensions, the preferred 
local ordering is icosahedral; this does not lead to a peri- 
odic crystalline structure and is the source of frustration. 
To model this type of frustration, we forbid certain local 
arrangements of the particles on the lattice. This can be 
done using two or n-body interactions. We follow [5] and 
set the interaction energy to be infinite if a particle has 
strictly more than £ particles as nearest neighbors. The 
interactions thus act as "geometric" constraints. 

Note that these geometric constraints are very similar 
to the kinetic constraints of the Kob- Andersen model [3]. 
However, as we shall show, the fact that they are encoded 
in an energy function makes a big difference physically, in 
particular for the thermodynamic behavior. 

We work in the grand canonical ensemble and intro- 
duce a chemical potential. All energies being zero or infin- 
ity, temperature plays no role. The thermodynamics for a 
given system (i.e., a given lattice) is then defined by the 
grand canonical partition function 

3(n) = ]T C e (n 1 ,...,n N )e^=^. (1) 

rai,...,nj\re{0,l} 

The dynamical variables are the site occupation values: 
n, = if site i is empty and rn = 1 if a particle is on that 
site. We take all the particles to be identical. In Eq. (1), ^ 
is the chemical potential, and CV(ni, . . . ,tin) implements 
all the geometrical constraints: it is the product of local 
constraints, one for each of the N sites. The term for site 
i is 

eii-m ]T nA (2) 

V jeV(i) / 

where Af(i) denotes the set of neighbors of i, and 8(x) = 
if x < 0, 9{x) = 1 if x > 0. 

For sake of concreteness and simplicity, we will focus in 
the core of the paper on the t = 1 case when each particle 
has at most one neighbor, deferring the general I case to 
Appendix A. Most numerical results will be given for the 
"basic model", noted BM, for which I = 1 and k = 2. 



2.2 Observables 

For a given "lattice" of N sites and type (Euclidean, tree, 
random graph, . . . ), and for a given form of constraints 
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(i.e., a value of £), there is only one parameter, the chemi- 
cal potential p. It is useful to introduce the grand potential 
and its density ui(p) = Q(p)/N, so that 

E{p) = exp [— pfi(p)] = cxp [— Npu>(p)] . (3) 

The pressure is given by p = —pui(p) and the particle 
density is 

/ x _ / E» n i \ = _J dSjp) = djpuj) = dp{p) 

PW ~\N/ NS(ji) dp dp dpi ' 

(4) 

Clearly, p is an increasing function of p. When p — > — oo, 
the system becomes empty a typical equilibrium configu- 
ration having almost all iii — 0, so p — > 0. In the opposite 
limit, /i — > +oo, there is in effect a strong penalty for each 
vacancy, but the geometrical constraints prevent all sites 
from being occupied; then p has a maximum value strictly 
less than 1. When N — > oo, this value is expected to con- 
verge to a limit p^. Obviously, p^ depends on the type 
of lattice and on the parameter I. 

We can define similarly the entropy density s(p): 

= ln - MaO = -^ip) - MaO (5) 

where again s(/z) should have a well defined thermody- 
namic limit. 

Other physical quantities that we shall study include 
susceptibilities associated with two-site connected corre- 
lation functions of the type (mnj) c - At low densities the 
m have only short range correlations. When p increases, 
p also increases; then the constraints become more impor- 
tant and correlations grow. When the density is close to 
Poo, the system will be very "rigid", allowing few fluctua- 
tions in the local density. It is plausible that the suscep- 
tibility will diverge at some critical value of p, separating 
a liquid phase at low p from a denser phase at large p. 

The nature of this dense phase, and of the transition, 
will depend on the lattice, on the boundary conditions, 
and need not be associated with a crystalline order. When 
it is not crystalline, we want a statistical description of 
the dominant equilibrium configurations. In particular, if 
these configurations form clusters, it is of interest to es- 
timate the number of such clusters. We will do this by 
computing the "configurational entropy" S(p) (also called 
complexity) associated with the number of clusters (also 
called states) of configurations with a given density p. 

3 Models on Cayley trees 

3.1 Iteration equations 

We first consider glass models defined on regular trees, i.e. 
connected graphs with no loops and fixed connectivity. 
We distinguish rooted trees where a site called the root 
is connected to only k neighbors while all the other sites 
(except for those at the surface, that is the leaves of the 
tree) have k + 1 neighbors. A Cayley tree is obtained by 
connecting a site to the roots of k + 1 rooted trees. 




i 

Fig. 1. An iterative method is used to compute the partition 
function on rooted trees. We begin with k rooted trees with 
roots j = 1, . . . , k and form a new rooted tree by joining a site 
i to each of them. The possible occupations of site i depend 
on the occupations of the sites j — 1, . . . , k and on the type of 
constraint used. 

When the graph is a rooted tree, the grand canonical 
partition function can be computed by recursion starting 
from the leaves. To do this, we follow conditional parti- 
tion functions because we need to know how to apply the 
constraints when joining the sub-trees (cf. Fig. 1). This 
provides a generalization of the well-known transfer ma- 
trix method for one dimensional systems (which can be 
viewed as rooted trees with k = 1). 

In our class of constraints, we need to know whether 
the root sites are occupied, and if they are, whether they 
have I occupied neighbors or less than that. Let ^\ , 

and be the conditional partition functions for 
a rooted sub-tree i when its root node i is empty (e), 
occupied but the constraint unsaturated (u) and finally 
occupied and the constraint saturated (s), i.e., the root 
site has I neighboring particles. Then the conditional par- 
tition functions for the rooted tree obtained by joining the 
sub-trees are easily computed. 

For instance, in the I = 1 case where each particle can 
have at most one neighbor, when merging k rooted trees 
(j = 1, ... k) to obtain a new tree rooted say at site i, we 
have 

-f^n-f, (7) 

Naturally, the total partition function is the sum of the 
conditional partition functions. 

To study these recursions, it is convenient to consider 
local fields defined via ratios of conditional partition func- 
tions. Here wc introduce on a root site i two fields a; and 
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hi defined as 



-(e) 



-fidi — 



~( e ) 



=■(«) i =•(«) ' 
"r" 1— 'i 



(9) 
(10) 



The first quantity is the probability for the root of a rooted 
tree to have an empty site (e); the second is the ratio of 
that probability and the probability to have a non satu- 
rated site. The use of /x when defining these fields is to 
simplify the analysis at large /x (cf. Sect. 5). These two 
fields form a closed recursion under the joining of sub- 
trees; for instance when 1 = 1, 



at = a(ai,bi, . . . ,a k ,b k ) 
= Iln 1 + e M(i-E 

bi = b(ai,h, . . . , Ofe, b k ) 
1 




1) 



= - In 



1 _|_ e f*(l-£;=i a*) 



(11) 

(12) 



Note that the use of ratios of conditional partition func- 
tions leads to recursions for two quantities rather than the 
initial three. To recover all the information in the initial 
recursions, we also keep track of the change in the grand 
potential. If . . . fi k give the grand potentials of the k 
sub-trees, we have after the merging 



+ Af2 itei (ai,bi, . . . , Ofe, bk) (13) 



where Afii tei (ai,bi, . . . , a k , bk) = Af2i tei is defined via 
E- E (e) + E {u) + E (s) 



-jj,AQ\i 



(14) 

With our definition of the fields, we have then ziAter = 
—a . From here on, we shall use the short-hand notation: 
h = (a, b) and hi = h(h\, . . . , h k ) with 



h (hi = (on, h), ...,h k = (a k ,b k )) = 
(a(ai,bi, . . . ,a k ,b k ),b(a\,b\, ■ ■ -,a k ,b k )j 



(15) 



3.2 Liquid phase 

Begin now on the leaves of a rooted tree, assuming that 
some kind of boundary condition is specified there. For 
example, the rij could be fixed, or their probability distri- 
bution could be given. These determine the initial values 
of the conditional partition functions and thus of the fields. 



We iterate the recursions, propagating the fields away 
from the leaves by performing mergings. When /j <C — 1, 
these iterations are contracting and so the fields converge 
to a value that is independent of the starting values on the 
leaves. The distribution of fields in the bulk (away from 
the leaves) is then given by 

V(h)=S(h-h Viq ) (16) 

with the fixed point condition h\i q = h(h\iq, . . . , /iu q ). We 
determine the fixed point for all /x, and refer to it as the 
liquid solution; its physical relevance includes at least the 
region \i <C — 1. 

We can now merge consistently k + 1 rooted trees to 
obtain a Cayley tree. We then have a liquid (or param- 
agnetic) phase, all correlations being short range and the 
heart of the Cayley tree being insensitive to the boundary 
conditions, even though a finite fraction of the sites live 
on the surface. In this regime, the homogeneous interior 
of the Cayley tree can be used to define the Bethe lattice 
model. 

In this context, also known as the Bethe- Peierls ap- 
proximation, the density u> can be obtained from the fol- 
lowing construction. Start with (k + 1) Bethe lattices; 
for each of them, pick an edge and remove it, leading to 
2(k + 1) infinite rooted trees. Then form two Bethe lat- 
tices by adding two sites and connecting each to (k + 1) 
of the rooted trees. Now the difference in grand potential 
between the resulting two Bethe lattices and the initial 
ones is just twice the grand potential per site (since two 
sites were added) and can be written as 



2w 



-(k + l)An edgc + 2A[2 sit 



(17) 



where Z\J? d g0 is the difference in grand potential corre- 
sponding to adding an edge and Af2 s i ie to merging (k + 1) 
branches into a new site. Such quantities are easily ex- 
pressed with the partition functions of rooted trees; for 
instance for I = 1 we obtain 

g-^^site(ai ,l>i>-",a*+i>6*+i) 

= i + gMd-sji 1 «,) L + j2 (e ^ _ 1} ^ (18) 



and similarly 



-llAOedge (01,61,02,62) 



-(e) -(e) -(e) /„(„) - 
— 1 —2 ' "1 I "2 



~(e) 

2 



"1 "2 



ai+a 2 ) ^gM a l 



(19) 



Note that we have the simple relation 
AQ s itc{h\, . . . , /ife+i) =AQ- lteT [h\, . . . , h k ) 

(h(hi,. 



, h k ), h k+ i^j 
(20) 
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Close to a liquid solution we have to first order 



Fig. 2. p dependence of the particle density p and entropy s 
of the liquid phase for the BM (£ = 1, k = 2). The failure of 
the liquid phase to correctly describe the model at high p is 
evidenced by the negative sign of the entropy for p > p s= o ~ 
7.40. But in fact the liquid solution becomes linearly unstable 
well before, at /i mo d — 2.77 (vertical line). 



whose interpretation is clear: the addition of one site by 
merging fc+1 rooted trees T\ , . . . , Tk+i can be decomposed 
into the construction of a new rooted tree T obtained by 
merging a site to the k rooted trees T\ , . . . , T k and the 
addition of an edge between Tq and Tk+\. Moreover Z\i? s ite 
is obtained from Z\Atcr by making the substitution k — ► 
k + 1. The expressions (18-19) are written for the general 
inhomogeneous case but simplify in the liquid phase where 
all the fields take their liquid value. 

In Fig. 2, we show as illustration the liquid's density 
p and its entropy density s, as a function of p for the 
BM. Since the models are discrete, the equilibrium en- 
tropy should never go negative. Nevertheless, the liquid 
solution at large p,, p > p s =o, leads to sn q < except in 
very special cases such as k = 1. Thus there must exist a 
phase transition at p c < pi s =o, and the liquid phase can- 
not be the equilibrium phase at /i > p c . Clearly, one must 
determine when the liquid solution is physically relevant; 
when it is not, one should find other solutions [5]. 



3.3 Linear stability limit of the liquid 

The Bethe-Peierls approximation fails when the bulk prop- 
erties of the Cayley tree become sensitive to the boundary 
conditions. This "instability" may show up via a loss of 
stability of the fixed point equations as given by a simple 
linear analysis. Indeed, starting with fields identically and 
independently distributed on the leaves, the field distri- 
bution V g +i(h) at "generation" g + 1 is related to that at 
generation g by 



/, re 
YldhjVgih^S (h i -h{h 1 ,...,h k ) 



i=i 



(Sh) g+1 = k^- 



(5h) t 



(22) 



liq 



where Sh = h — hn q and (.} g refers to the average using 
the distribution V g . Since ft is a two-component vector, 
dh/dhi is actually a 2 x 2 Jacobian matrix. If Ai denotes 
the eigenvalue of largest modulus of that matrix, the sta- 
bility criterion simply reads 



fe|Ai| < 1. 



(23) 



When fc|Ai| > 1, the liquid solution is unstable to pertur- 
bations homogeneous within a generation; we shall refer to 
it as the modulation instability because it is a transition to 
a regime with successive (homogeneous) generations car- 
rying different fields. 

An alternative point of view consists in studying re- 
sponse functions. At finite p, the response to a pertur- 
bation is related to correlations through the fluctuation- 
dissipation theorem, and an instability can be detected by 
means of a diverging susceptibility. Thus, we recover the 
previous result by computing the linear susceptibility in 
the liquid phase, looking for the point where it becomes 
infinite. The linear susceptibility is defined in terms of 
connected correlation functions {niUj) c via 



(24) 



Making use of the homogeneity of the liquid solution, it 
can be rewritten 

oo 

Xi (/i) = P(l -P) + J2( k + WH^n r ) c (25) 

r=l 

where n and n r are taken at distance r in the tree. The 
series converges provided that 



mfc+ l im ln <"°^ < Q. 



(26) 



To evaluate (non r ) c , we invoke the fluctuation-dissipation 
relation 

(n a n r ) c - (27) 



(c) 

where h denotes the external field conjugate to no- Since 

/iq C ' is a function of (the components of) h only, we can 
use the chain rule 



fr—\ 



d(n r ) d(n r ) ( j-j dhi \ dh 



dh ( c} dhr-i {fJi 9^-! J dh ^ 



(28) 



(21) 



where we introduced hi = h(hi-i, hn q . . . , /ii; q ) as inter- 
mediate fields. In the liquid phase, all these fields are 
equal and the previous equation factorizes, leading again 
to k\Xi\ < 1. For instance for the BM, the modulation 
instability shows up at ^. mo d = 4 In 2 ~ 2.77, well before 
the entropy becomes negative at /i s= o — 7.40 as shown in 
Fig. 2. 



6 



O. Rivoire et al.: Glass models on Bethe lattices 



3.4 Crystal phase 

We can ask whether it is possible to choose boundary con- 
ditions such that the interior of the Cayley tree has a peri- 
odic structure (for a general tree, the boundary conditions 
will vary from leaf to leaf). In the /x — > oo limit, we expect 
a crystalline phase to exist whose structure can be easily 
displayed by starting at the center of the Cayley tree, fill- 
ing the sites with particles whenever possible. In this case, 
for the BM, the field h = (a, b) takes three different values 
h e , h u and h s such that 

h e = h(h s ,h s ), 

h u = h(h e , h e ), (29) 
h s = h(h e ,h u ), 

with h e = (0,0), h u = (1, 1) and h s = (1,0). The integer 
values and 1 of the components of these fields reflect the 
fact that a given site is certainly empty, unsaturated or 
saturated. For finite fi, fluctuations are present but we can 
still look for a fixed point of Eq. (29), with h e ,h u ,h s £ R. 

Such a solution, distinct from the liquid one, is found 
to exist and to be stable for \i > fi ms with p ms ~ 2.89 (ms 
for melting spinodat). Next we want to evaluate the grand 
potential of this crystalline solution in the bulk. To do so, 
we first consider the \i = oo limit and estimate the density 
of empty and occupied sites, po = 2/5 and p\ = 3/5, and 
also the proportion of edges connected to one and two 
particles, tti — 4/5 and 7T2 = 1/5. Then the crystalline 
potential can be written as w crys t = Lo S i tc — 3/2 Wedge, with 

t^sitc = poAf2 aite (h s ,h s ,h s ) + piAfi site (h u ,h e , h e ), 

Wedge = 7Tl^^edge(fre, h s ) + TT2 A f2 c d gc (h u , h u ). 

Now for p ms < [i < oo, the structure of the crystal is 
preserved, with empty and occupied sites being replaced 
by most probably empty or occupied site. Therefore we 
resort to Eq. (30), using the adequate values of the fields 
h e , h u , h u , and we keep the same factors po,i an d 7To,i. 
This leads to a melting transition (where wii q = u cryst ) at 
Mm 3.24. 

Comparing fi ms ~ 2.89 with the location of the liquid's 
instability p mo d — 2.77, we have an interval ^ mo d < t 1 < 
[i ms where no homogeneous nor periodic solution seems 
to exist. This is to be contrasted with the phase diagram 
of the hard sphere model studied by Runnels [17] where 
the presence of a particle on a site forbids the occupation 
of any of its neighboring sites. That model has been re- 
considered recently in two dimensions as a combinatorial 
problem of counting binary matrices with no two adjacent 
l's when \i = [18] and on random graphs as an optimiza- 
tion problem called vertex cover problem when fi = oo [19]; 
from the point of view of our lattice glasses, these models 
correspond to the £ = constraint. In this case the modu- 
lation instability coincides exactly with the melting tran- 
sition (and with the spinodal point), fi m = fj, ms = /i mo d- 
This is due to the special structure of the crystal, or- 
ganized in alternate shells of empty and occupied sites, 
and therefore described by a cyclic solution of the liquid 



equation, h\ = h(h , . . . , ho) and ho = h(h\, . . . , h\) with 
homogeneous shells (i.e., the arguments of h are all iden- 
tical). For £ = 1 no such cyclic solution was found and 
homogeneous boundary conditions on the leaves yield an 
aperiodic behavior of the recursion hj + i = h(hj, . . . ,hj). 
This feature is very specific to the unphysical nature of 
pure Cayley trees: it does not survive in the random reg- 
ular graphs which we use below. Therefore we have not 
pushed its study any further. 

4 Models on random regular graphs 

4.1 From Cayley trees to random regular graphs 

When the Bethe-Peierls approximation no longer holds on 
a Cayley tree, the sensitivity to boundary conditions does 
not allow one to define a thermodynamic limit and there- 
fore may lead to unphysical results. Since that is due to 
the presence of a finite fraction of the sites on the surface, 
one way to get rid of this problem is to define Bethe lat- 
tice models as models on random regular graphs. These are 
simple graphs with fixed connectivity k + 1, simple mean- 
ing that they have no trivial loops (joining a site to itself) 
and no multi-edges (no two edges join the same sites). 
Here we will use the cavity method to obtain results for a 
typical random graph (chosen uniformly in the set of all 
random regular graphs with a fixed connectivity k + 1). 

When the number N of sites is large, typical random 
regular graphs look locally like trees, having only long 
loops of order In N. Therefore the recursive equations still 
(locally) hold, making these lattices analytically tractable. 
The large loops implement an analog of generic boundary 
conditions and the resulting frustration forbids crystalline 
orderings. One then expects the system to possess an equi- 
librium glass phase in the high /i region. 

In the low /i phase, the liquid solution is recovered. The 
corresponding Bethc-Pcierls approximation is called the 
factorized replica symmetric approximation in the context 
of the cavity method where the vocabulary is inherited 
from the treatment of spin glasses based on the replica 
trick. However the glassy high \x regime will be character- 
ized by the existence of many solutions of the local equa- 
tions and we will have to resort to the replica symmetry 
breaking (rsb) formalism to correctly take into account 
the specific organization of these solutions [20]. 

4.2 Entropy crisis 

Reconsidering the stability of the liquid with the modu- 
lation instability now being excluded, we can look for a 
"spin-glass" instability, to borrow the terminology from 
magnetic systems where the modulation instability is re- 
ferred to as the "ferromagnetic instability" . This new in- 
stability manifests itself as a divergence of the non-linear 
susceptibility, which is defined as 

*2M = ^E^)c- ( 31 ) 

i,3 



O. Rivoire et al.: Glass models on Bethe lattices 



7 



In the formalism of Sect. 3.3, the instability appears as a 
widening of the variance ((Sh) 2 ) g under the recursion of 
Eq. (21). Both approaches are equivalent and lead to a 
stability criterion 

fc|Ai| 2 <l. (32) 

The eigenvalue Ai is the same as for the linear suscepti- 
bility since the transfer matrix is just the square of the 
Jacobian matrix we had in Eq. (22). Note that this con- 
dition is always weaker than that for the modulation in- 
stability fc|Ai| < 1. However it is the relevant one in the 
case of random graphs where homogeneous perturbations 
are generically incompatible with frustration. If the liq- 
uid is locally stable for all /z, a continuous phase tran- 
sition is excluded. For i = 1 this happens for k = 2, 3 
because Vk\\i(fj,)\ < 1 for all fi with only asymptotically 
\/fc|Ai(/i)| — > 1 as (i — > oo. In the general case, as soon as 
[i g > Ms=0: where ji g (g for glass) is defined by 

\/fc|Ai(M = M 3 )| = 1, (33) 

the resolution of the entropy crisis requires a phase tran- 
sition before the spin-glass local instability is reached, 
and we conclude that a discontinuous phase transition 
must take place at fi c < /i s =o- I n that case we expect 
a behavior similar to that of infinite-connectivity models 
solved within a one-step replica symmetry breaking (1- 
rsb) Ansatz, like e.g. the p-spin models (p > 2). When 
p g < Ms=Oi as is found for I = 1 and k > 4, we can ei- 
ther have a continuous transition at /i g or a discontinuous 
transition at /i c < ji g ; a study of the local stability of the 
liquid solution says nothing about which case arises. 

4.3 Cavity equations 

The solution by the cavity method predicts results for 
quantities averaged over all random regular graphs with 
size N — > oo, but the problem is more clearly stated on a 
given finite regular graph. Indeed, we want to solve self- 
consistently a set of (k + 1)N coupled equations for the 
cavity fields 

hi^jo = M^ji— >*i • • • j hjk^i) (34) 

v« v{io,ii,.-.,ife} e Af(i) 

where Af(i) denotes the set of k+1 neighbors of i. Thus for 
each i £ {1, . . . , N} we have k + 1 equations correspond- 
ing to the different choices of jo € Af(i). The notation 
hi^j refers to the local field on i when the edge to the site 
j is removed, and it is called a cavity field; for a Cayley 
tree, it corresponds to the field on the root i of a rooted 
tree obtained when the edge ij has been removed. These 
local equations are known as the TAP equations in statis- 
tical physics [14] and as the belief propagation equations 
in computer science [21]. The contact between both points 
of view has recently lead to establish that these equations 
must always have at least one fixed point corresponding to 
the minimum of a correctly defined Bcthe-Peierls approx- 
imated free energy [22] . Their solutions should correspond 



to the fuzzy concept of states ubiquitous in the spin-glass 
literature. 

ilgorithm can be used to try to 
solve these equations on a given graph of large but fi- 
nite size N. First on each oriented edge we associate local 
fields hi^j randomly initialized. Then we proceed itera- 
tively: at each step, all the oriented edges are successively 
chosen in a random order, and the field on the chosen ori- 
ented edge is updated by taking into account the values 
of its neighbors as prescribed by Eq. (34). The iteration 
is stopped when a sweep of all the oriented edges results 
in no change; in such a case we lie at a fixed point of 
the local equations. However the algorithm may also not 
converge, in which case no conclusion can be drawn. In 
practice, we find a rapid convergence toward the liquid 
solution for p < m, p and then a failure to converge for 
A* > Mftp- The critical /ibp depends slightly on the graph 
but for large N and large k, it is given by the glass in- 
stability, i.e., [itp — *• Ms as k — > oo. In some case (for not 
too big graphs), we could find a convergence towards a 
non-liquid distribution, suggesting that the high /i region 
corresponds in fact to a glassy phase. In order to deal with 
this phase where the Bethe-Peierls approximation breaks 
down and simple message passing algorithms fail, we will 
now introduce the cavity method which provides both an 
alternative approximation for infinite graphs and insights 
into elaborating more efficient algorithms for finite graphs 
[16]. In Appendix D we present an alternative approach 
based on the replica method. 

4.4 One-step replica symmetry breaking cavity method 

The local equations Eq. (34) can be written for arbitrary 
graphs but they provide exact marginal probability distri- 
bution (and thus exact particle densities) only for trees. 
In addition, they are in general intractable. However, they 
are particularly suited for very large random graphs, where 
due to the local tree-like structure, they are expected to 
provide good approximations of the marginals and where 
additional hypotheses allow for an analytical treatment. 
To do so, we do not try to find one solution but instead 
turn to a statistical treatment of sets of solutions. 

Being interested in the case where many solutions ex- 
ist, we fix a /i for which this is supposed to happen. We 
make the further hypothesis that exponentially many (in 
N) solutions exist. More precisely, we assume that the 
number A/at(w) of solutions with a given potential density 
u> on graphs of size N is given by 

Af N (u) ~ exp[NS(uj)} (35) 

where S(u) > is called the configurational entropy (or 
complexity) and is supposed to be an increasing and con- 
cave function of the grand potential w. This is a strong 
hypothesis which is justified by its self-consistency and by 
its consequences (in particular it matches with the output 
of replica theory calculations, cf. Appendix D). 

Starting with a graph Q, we pick a site i and one of 
its neighbors jo <= and define the graph Gi^j as 
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the connected graph containing i obtained by removing 
the edge ijo from Q. If Q is a Cayley tree, Gi^j is nothing 
but a rooted tree, the appropriate structure to write down 
recursive relations. We introduce Ri—>j (h,u>), the joint 
probability density, when a solution a of Eq. (34) defined 

on Gi^j a is chosen randomly, that the cavity fields 

on i take the value h and the grand potential density - o 
the value w, 

«) = ^ £ 5 - h ^ ) s (- - -So) • 

b a 

(36) 

The next crucial hypothesis is based on the locally tree- 
like structure of large random graphs. Indeed, if the graph 
Gi^j is a rooted tree, the sub-rooted trees Gj 1 ^i 1 ■ • • , Gj k ^i 
with {ji, . . . , jk} € A/"(i) are disjoint and the cavity fields 
hj 1 ^i, . . . , hj k ^i are therefore solutions of uncoupled local 
equations. In such a case the joint probability distribution 
factorizes over the independent variables 



consider only potentials to close to u>q, noted lu e V Wo , 
such that we can linearize the complexity 



with 



Jl—nV<-Jl—ni ^Ji^r 



,R 



(37) 



In addition, a simple relation between the distribution 



on 



Ri^jo on a rooted tree Gi^j and the Rj 1 ^i,. . . ,Rj k - 
its sub-rooted trees Gj^i,- ■ ■ ,Gj k ^i can be written: 

Ri^ jo (h,uj) = / dh jr <L) jr R jr ^i(h jr ,u) jr ) 

■* r=l 

5 (ft - h({ V») ^ W - 1 ^ JV> jp + Mter({/l jV })^ 

(38)' 

where the A^ r are the sizes of the sub-trees and N = 
J2 r Nj r + 1 is the size of Gi^j ■ A further simplification 
takes place on regular trees where, due to the absence of 
local quenched disorder, the equations for hj 1 —^i,. . . , hj k ^i 
are identical. In such a case, Ri^j is in fact independent 
of the oriented edge i — ► jo- For a large random regular 
graph, we assume that the same properties hold. Note that 
the last hypothesis, yielding Ri^j a = R is called the fac- 
torization approximation and could be relaxed by working 
with a distribution TZ[R] of the R over the various edges 
[20]. However the factorization approximation should be 
exact on random regular graph for systems without dis- 
order provided that no spontaneous breaking of "trans- 
lational" invariance occurs, and we will not consider this 
extension here. 

We now want to write the 1-rsb cavity equation, which 
is a self-consistent equation for the distribution P^ of lo- 
cal fields h at a fixed grand potential u. P^ is obviously 
proportional to R, P^\h) oc R{h, lu), and since the distri- 
bution of the to is given by / dhR(h, lu) = C exp[iV S '(ui)] , 
we have the relation 



R(h,w) = Ce* TS MpM(h) 



(39) 



with C a proportionality constant independent of both 
h and to. Next we fix a grand potential density lu and 



m(w ) = -zr-{uo)- 
H duj 



(40) 



(41) 



Given the concavity of the complexity S, fixing luq is 
equivalent to fixing m and thus V UQ can be rewritten as 
V m - Plugging relation (39) into Eq. (38), we find that the 
distribution defined by 

P (m) (/j)oc f duuP { ^(h) (42) 

satisfies a simple self-consistent equation 
P {m) {h) oc 

/ YldhjP^ih^S (h-h({hj}j) e -™^({ fc 3 }). 

(43) 

Eq. (43) is called the factorized 1-rsb cavity equation. We 
will drop the explicit reference to the parameter m in the 
ensuing discussion, but it should be kept in mind that the 
1-rsb cavity field distribution P(h) is m-dependent. 

As we have seen, it is convenient to fix m instead of u>; 
going from lu to m actually amount to performing a Leg- 
endre transformation. The complexity S(co) is recovered 
by Legendre transforming the function <P(m) defined as 



m<£(ra) = muu 



via the relation 



1 

H 



S{uj) = m 2 d m <P{m). 



(44) 



(45) 



Eq. (44) is similar to the definition of the entropy 
through 

fiuj(fi) = ~np- s(fj,) (46) 

[cf. Eq. (5)]. Indeed <P(m) is for the states the analog of the 
grand potential to{fi) for the configurations, the parameter 
m sampling the states as the chemical potential /j, samples 
configurations, and the complexity S(cu) counting states 
as the entropy s(fi) counts configurations [23]. The quan- 
tity corresponding to the grand partition function is 

E(m) = ^ exp[-iVm^ (a V)] 

= J duj cxp (N[E(u) - mfiLu}) ( 4? ) 

= exp[— Nm^{m)}. 

Following the same analogy, <P(m) can be computed simi- 
larly to a grand potential. To do so, we first need to gen- 
eralize to random regular graphs the construction that led 
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us to the Bethe-Peierls approximation of the grand poten- 
tial, Eq. (17). Likewise, we want to write 



<2>(m) = A$ Bite (m) 



1 



with A(I> s it e (m) the contribution from a site addition, and 
A<P c d ge (rn) from an edge addition. The way two sites can 
be added is even simpler than for Bethe lattices here be- 
cause we do not care about introducing loops: take a ran- 
dom regular graph of size N and connectivity k + 1, pick 
(k + 1) edges and remove them, leading to 2(fc + 1) ampu- 
tated sites with k neighbors instead of (k + 1). Then add 
two new sites and connect each one to (fc + 1) of the ampu- 
tated sites, leading to a new random regular graph of size 
7V+2 and same connectivity (fc + 1). Thus the contribution 
to <P when going from N to N + 2 sites is equivalent to 
the contribution from two site additions plus (fc + 1) edge 
deletions [i.e., minus (fc + 1) edge additions], as expressed 
by Eq. (48). 

To see how this is related to the grand potential shifts 
Z\J? s itc and Z\j? dg , we rewrite Eq. (39) introducing the 
definition Eq. (41) of the parameter m as 



R(h,u) =e m ^ N ^-^ 0) )p{h) 



(49) 



where uj^ = u)(°\m) enforces the normalization. The 
function 



O m^(f2-f2 (0) ) 



gives the distribution of grand potential [2 = Nlu on a 
rooted graph of size N with one site having only k neigh- 
bors. Now we take fc+1 such rooted graphs and determine 
the distribution p\{Q) when one site is added. It is given 
by 



Pl {(2) 



fc+1 

\ \ dhjdOjRihj^^S 
3=1 



= e 



fc+l 



(o)\ 




We compute similarly the contribution from edge addi- 
tion, 



APedgeM (48) P2(rt) 



f 2 I 2 

= / dhjdf2jR(hj, f2j)8 I fl — flj — Af2 c a gc (hi, hz) 

2 

= e M n - n ™-°™) [ f[dh j P(h j )e- m ^^ h ^ 

J 3=1 

_ e m^(j"2-r2< 2 )) 



(53) 

with {2^ = Q[ 0) + l^ 0) + A$ cdgc , while the mean shift 
^^odge due to an edge addition is 



I r 

A<P edgc (m) = In / Yl dh j P(h j )e- m " Aned ^ hl ' h ^ 

nifj, 

(54) 

The replica symmetric description of the liquid phase 
is recovered by taking P(h) = 6(h — /in q ). When many 
solutions coexist, by varying m = m(ujo) at fixed ^, we 
describe states characterized by different values of u)o and 
the question is which m must be selected to describe the 
equilibrium (glassy) thermodynamics. Following Eq. (47), 
the grand partition function is 



= J dujcxp (N[S(w) - fiu]) 



(55) 



with = rtj 



and the saddle point method for N — > oo indicates that 

(50) 

the grand potential u> of the dominating states is such 
that jU = d w E(w). But of course, this saddle equation 
is relevant only if its solution u s lies inside the interval 
range ]w m i n , w max [ where the integral is performed, which 
corresponds to the range where S(u>) is strictly positive. 
Generically, for \i > /j, s , it is found that uo s < iv m i n , and 
we must then take instead uj s = w m i n . A kind of replica 
trick intervenes here to balance the complexity contribu- 
tion thanks to the parameter m, in such a way that iv s 
is always given by a saddle equation. So, in the glassy 
^Asitc({%$)aKc fi > fi s where the rsb formalism becomes necessary, 
we Avant the equilibrium grand potential to be given by 
lo s — cj m i n (^) = <P(m s ). Since the complexity curve is ex- 
e({hj}) pected to be continuous at w m i n , we can alternatively ask 
for the condition £(ui s ) = 0. From Eq. (45), we see that 
it corresponds to extremizing <P(m), i.e., 



A$ s 



(51) 



while the mean shift 



d m <P(m 



,) = 



(56) 



Z\^ s ito due to a site addition is 



A$ site (m) = In 



fc+i 

YldhjPih^e-™^* 

3 = 1 



e({hj}) 



(52) 



which is precisely the criterion provided by the replica 
method for selecting the breaking point parameter m. 

Note that other values of m also carry physical infor- 
mation. Lower values of m (m < m s ) describe metastable 
states for which w > w s (it is the analog of a non-zero tem- 
perature giving access to excited configurations); of partic- 
ular interest is the value LOd associated to the maximum 
complexity S(u>d) = max w S(lu) since its describes the 
most numerous states. We expect that this is the portion 
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of phase space where the system will get almost trapped 
at long times after a "quench" from the low density liquid 
phase. 

Higher values of to (to > m s ) are usually termed as 
unphysical, but in fact they describe properties of systems 
associated with untypical graphs. In particular, it can be 
shown that to = 1 always gives back the liquid solution. 
In fact, for to > m s , the cavity method leads to a nega- 
tive complexity, which seems to be in contradiction with 
its initial definition Eq. (41). The point is that in the cav- 
ity method, the graph is not specified and e NE<<0J ^ is the 
number of states with grand potential u> after averaging 
over different graph realizations. This fact has no conse- 
quence whenever the average corresponds to the typical 
case, which is expected as soon as S(lu) > ; indeed some 
graph realizations may behave very differently from typ- 
ical realizations, but their contribution to the averaged 
quantities is negligible. However, an exception is worth 
mentioning: if the quantity we average is typically strictly 
zero, but happen to be positive for exponentially rare real- 
izations, it leads to a small but non zero (i.e. non typical) 
average. We find such a behavior here, where some untypi- 

,(typ) 



cal graphs allow for uj s lower than the typical value ojI 

leading to a S even if the complexity on 

a given graph is intrinsically a positive quantity. 

An analog phenomenon happens for instance in the 
random energy model, where averaging over disorder leads 
to a negative entropy associated with energies lower than 
the typical ground state, while for a given realization of 
the disorder the entropy is necessarily positive. Here the 
role of the quenched disorder is taken by the topological 
disorder from the various realizations of random regular 
graphs. From this point of view, a random graph with no 
frustrating loop is an example of untypical graph which 
has a crystalline ground state, as opposed to the glassy 
ground states of typical random regular graphs. 



4.5 Observables 



To complete our overview of the cavity method, we now 
show on the example of the particle density how physical 
observables can be computed; Sect. 4.9 will provide an 
other example with the computation of susceptibilities. 
We begin by considering rooted trees where the particle 
density on the root i is simply given by 



s (s) + s {u) 

(™i)rootcd tree = „( s ) ]_ „(„) * „( e ) = 1 ~ e ^ ■ (57) 
1 — ' i ' 1 — i ' ' — ' i. 



Now for a Cayley tree, we need to take into account k + 1 
neighbors instead of k, 



where the total local field A is computed similarly to the 
cavity field a but substituting k by k + 1, i.e., 



Ai = A(ai,bi, . . . ,afe+i,6fc+i) 



fe+i 



l + e^i-E*;^) (l + ]T(e^-l) 

i=i 



(59) 



The total field Bi would be defined similarly; note that 
with our notations, we simply have A — — AO s i te . On a 
random graph in the 1-rsb phase, we need to include the 
reweighting associated with the addition of the site i: 

p{H, to) 

JUjtl dhjP{hj)(l - e -M({M)) e -™MfW{M) 



/ Ujtl dhjP{hj)e- m » Afl ^ h i» 



(60) 



and the equilibrium value is given by p(/i, m s ). The same 
lines can be followed to compute any other observables. 



4.6 Order parameters 

Comparing with the replica method [14], the cavity ap- 
proach focuses on local fields instead of overlaps between 
states. In lattice glass models, the latter can be defined as 



1 N 

Qap-y... = y^ J {ni) a {n i )i3{ni) 1 . . . 



(61) 



with the indices a, (3, 7, . . . denoting states randomly cho- 
sen according to their Boltzmann weights. 

Overlaps are particularly useful in infinite connectiv- 
ity systems where, due to the central limit theorem, the 
first two moments q a and q a are sufficient to encode the 
Gaussian distribution of the cavity fields. In contrast, for 
finite connectivity systems, an infinite number of overlaps 
must be kept and working directly with the local field 
distribution is simpler. However, the two choices of order 
parameters, local cavity fields or global overlaps, provide 
complete and equivalent descriptions; in particular over- 
laps can be easily recovered from the knowledge of the 
cavity field distribution. 

At the replica symmetric level, all indices a, /?, 7, 
. . . are equivalent, so we only need to distinguish overlaps 
according to the number r of states they involve, and we 
have 



J dAdBV rs {A, B) (l - e-^Y (62) 



with A and B being the total local fields of Sect. 4.5. 
Note that the replica symmetric approximation in princi- 
ple already involves a functional order parameter V rs , but 
in our case where the factorization Ansatz is taken, the 



(^i)caylcy tree — 1 C ^ l . 



liquid distribution is trivial and we merely have qr 



(rs) 



(58) (1 — e Miiq)»^ j e j-jjg orc j er parameter is a single scalar. 
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At the one-step level, we need to distinguish whether 
randomly chosen states are distinct or identical. Thus for 
overlaps involving two replicas, we can define two param- 
eters, qo = q a f3 (a 7^ P) and q\ = q aai corresponding re- 
spectively in spin glasses to the Edwards- Anderson order 
parameter and the (squared) magnetization. Their com- 
putation amounts to calculating densities inside a state, 
which was the subject of the previous section. Note how- 
ever that go and 9i provide only a partial description of the 
system and the full order parameter has here a functional 
structure P(h). Any new level of replica symmetry break- 
ing will require considering a more sophisticated order pa- 
rameter, namely a distribution over the order parameters 
from the previous level. For instance, the two-step order 
parameter will be written as a distribution Q[P] over dis- 
tributions P(h). Describing with this formalism a finite 
connectivity system with full replica symmetry breaking 
is therefore rather complicated, and we will limit ourself 
to at most two levels of replica symmetry breaking in the 
ensuing discussion. 

4.7 Solution via population dynamics 

Given the cavity equations (43), we would like to solve 
them. Since they are essentially functional relations, an 
analytical treatment is not possible in general. An impor- 
tant exception however is the close-packed limit fi — > oo ; 
the next section is devoted to this case. Here we consider 
the more general finite /i situation and use the popula- 
tion dynamics algorithm of Ref. [20] to obtain numeri- 
cal results. The principle of the algorithm is elementary: 
the distribution P(h) is encoded in a family of M fields 
{/ij}i=i,...,M such that 

M 

^)^E^-U (63) 

i=l 

and the cavity field distribution is expressed as the fixed 
point of the iteration equation 

P g+1 (h)<x f Y[dh j P g (h j )5(h-h(h u ...,h k )J 

^ j=i (64) 

e -mnAn lter (hi,...,hk) 

At each step, M new children are generated; each is ob- 
tained by choosing randomly k parents among the pop- 
ulation. To take into account the reweighting, we dupli- 
cate or eliminate the children according to their weight 
e -mnAn itet SQ we keep a total population of (approx- 
imately) M individuals. 

However, such a recursion turns out to be unstable 
and we stabilize it by means of a relaxation parameter e e 
]0, 1[. At each step, only a fraction eM of the population is 
regenerated. The reason for this relaxation will be made 
clearer when we will discuss the stability of the cavity 
method solution; it will be associated to the instability of 
the first kind discussed in Sect. 5.2.2. 



4.8 Static and dynamical transitions 

At low n, the population dynamics algorithm always con- 
verges to the liquid solution, i.e., starting from a popula- 
tion with an arbitrary distribution it converges to a pop- 
ulation of identical fields corresponding to the fixed point 
/iiiq of h. When /x is increased, a first non-trivial distribu- 
tion is found at /i = fid for m = 1. At this point many 
states exist, but they are only metastable and the statics 
is still given by the (paramagnetic) liquid state. This er- 
godicity breaking is called a dynamical phase transition 
because it is where the equilibrium dynamics should dis- 
play an ergodic-non ergodic transition (see however the 
discussion in Sect. 6). 

A static phase transition, which is the one relevant at 
equilibrium, only appears for higher chemical potential, 
/i = ji s , when the configurational entropy vanishes. At this 
point the grand potential of the 1-rsb solution becomes 
lower than the one of the liquid and the equilibrium phase 
transition takes place. 

In practice, the dynamical transition point /id is found 
by decreasing /j, and looking for the \x where the 1-rsb so- 
lution disappears. For the BM, it happens at /id — 6.4 ; 
more generally for I = 1 and other k we find /id < Ms i indi- 
cating that the transition is always discontinuous, as could 
not be directly inferred from the arguments of Sect. 4.2. 
To obtain fi s , we calculate explicitly d m (j>(m, /i) and look 
for the /j, at which d m <f>(m = 1,/z) = 0; for the BM, we 
obtain n s ~ 7. 

4.9 Stability of the one-step solution 

To determine whether the equilibrium state is really de- 
scribed by a 1-rsb solution or whether further replica break- 
ings are necessary, one has to study the stability of the 1- 
rsb solution. In this section we set up the formalism needed 
to check it. In the following sections we will analyze the 
stability in the close packing limit (fi — > oo). 

In the 1-rsb phase the Gibbs measure is decomposed 
in a cluster of different thermodynamic pure states [14]. 
Thus, there are two different types of instabilities that can 
show up [14]. First kind: The states can aggregate into 
different clusters (see Fig. 3). To study this instability one 
has to compute inter-state susceptibilities: 

4 ntor = ^ E (WW -W)W)) P - (65) 

where the overlinc denotes an average over the states taken 
with their Boltzmann weights. Second kind: Each state 
can fragment in different states (see Fig. 3). To study this 
instability one has to compute intra-state susceptibilities: 

If any of the intra or inter-state susceptibilities diverge 
then the 1-rsb glass phase is unstable (toward a 2-rsb glass 
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Fig. 3. Pictorial view of the two possible instabilities of the 1- 
rsb Ansatz. At the bottom, the states stay states but clusterize 
(first kind) . At the top, the states become cluster of new states 
(second kind). If we were on a totally connected graph where 
an overlap function q(x) can be defined, its 1-rsb shape would 
be affected on different parts depending on which instability is 
relevant; indeed its left part (a; < m) corresponds to the inter- 
state overlap go and its right part (x > m) to the intra-state 
overlap q\. 




Fig. 4. Cavity diagram for computing a two-site correla- 
tion function {non r )c- As all we know is P(h), the distribu- 
tion of the local field on the root of a rooted tree, we build 
a chain of sites / = 0, . . . , r out of rooted trees with fields 
ho,gij, ■ ■ ■ ,g r -i,j,g r - We proceed recursively: we first add site 
1 = 1, obtain a new rooted tree with root 1 and local field 
hi = h(ho,gi t i, . . . , <?i,fc-i) so the length of the chain to com- 
pute is reduced from r + 1 to r. Then we proceed further by 
adding site 2, etc (see also Fig. 5). 



and 



((n ) (n r ) - p 2 ) - exp(-r/C 2 ). 



(69) 



Due to the tree-like structure of the lattice, the result- 
ing stability conditions are different from the condition 
£ < oo used for finite dimension lattices and reads 



phase). However, as for the liquid, the linear susceptibil- 
ities Xi are related to instabilities incompatible with the 
underlying random graph structure and, hence, they are 
irrelevant for our purposes. In the following we will focus 
on the p = 2 case which is the only relevant one since 
all the susceptibilities with p > 2 are clearly bounded in 
modulus by the p = 2 one. 

Because of the homogeneity of the simple random graphs 
that we are focusing on, the stability analysis is simplified 
and in particular: 



6< hTfc 



C2< InT 



(70) 



In the following we show how the correlation length £ p 
can be computed. A very similar procedure can be carried 
out for £2- In Appendix B we shall show explicitly how 
this can be done in the close-packing limit (fj, = 00). 

In order to obtain £ p we need to compute the corre- 
lation functions (no« r )?. We write the generalization of 
Eq. (28) as 



x^ ntor = E( fc + i)^ 1 ( - t 

oc 



(n a n r ) c 



(67) 




n °^{h, 



dhn 



dh, 



(c) 



(71) 



r=l 



where n and n r are at distance r (we omitted the unim- 
portant r — term). 

We expect, as it can be proved (see below), that the 
correlation functions decay exponentially at large distance 



(n n r ) 2 ~ exp(-r/£ 2 ), 



(68) 



where the fields arc (d — l)-dimensional vectors (d = 3 
for our models), so that the product involves in fact (d — 
1) x (d— 1) dimensional Jacobian matrices d\h (the nota- 
tion d\h indicates that the derivative is taken with respect 
to the first field, here and ()* means transposed). 

Note that h(hi—i;gij) is used as a short-hand notation 

for h(hi-i, gi.i, ■ ■ ■ ,gi,k-i) (see Fig. 4). Next we need to 
take into account the reweighting introduced by the addi- 
tion of the sites I = 1, . . . , r — 1. With a transfer-matrix 
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Y r -'> T7M 

Fig. 5. Partition functions involved in the reweighting for 
the comp utation at the 1-rsb level of the two-site correlation 
(n n r )S, c , illustrated here with k = 2. We start with the parti- 
tion functions S w ,T W , . . . , T^ 1 ', E (r) where T (l) is in gen- 
eral the product of the partition functions of k — 1 rooted trees 
(top of the figure); for k — 2 as illustrated, it reduce to one 
rooted tree with root site noted I' . Next we merge the k rooted 
trees corresponding to 5^ and T^ 1 ' into a site 1 and call the 
partition function of the resulting rooted tree (bottom). 
Recursively, we define similarly the other S"' for 2 < I < r — 1. 



approach in mind, we write it as 



~(o) (nr=i r(,) ) ^ (r) 



r-l 



1=1 v ' 



(72) 



where the notation refers to Fig. 5: denotes the prod- 
uct of the k — 1 partition functions associated with the 
rooted trees with cavity fields gij (1 < j < k — 1) and 
(1 < / < k — 1) the partition function of the rooted 
tree obtained by connecting a new site I to the k rooted 
trees corresponding to and 

We can write for one site addition 



~(0 



exp[-m/iZiAter(/i*-i; 9i,j)}, (73) 



and for the last edge addition 



cxp[-m/j,An cdgc (h r - 1 ,g r )}. (74) 



The 1-rsb formula for {non r )% is therefore 

>--i / fe-i 



Z- 1 J dh P(h ) 

d(nr)Y ( T ^ dl f l{hl _ l . gi 



i=i \ji=i 

dh 



dg r P{g r ) 



sl = l 



dh r -i / 

r-l \ 



(75) 



where the normalization Z is given by 

r-l / fc-1 



: J dh P(h 



dg r P{g r ) 



/r-l 



s/=l 



(76) 



To be complete, we also need to insert in the previous 
formulae the following identity defining the intermediate 
fields hi, 



1 = 



.r-l 

/ Hdh^fht -h(hi-! ;g hj )). (77) 



To determine the behavior of the correlation functions 
between sites at distance r, we introduce two transfer ma- 
trices, corresponding respectively to the numerator and 
denominator of Eq. (75), 



. fe-i 



T n {hi-x,hi)= / II dgi^PimMhihi^ 

S (h t - hiht^gij,)) e - m » A ^ h >-^i\ 
r fe-1 

T d (hi- U hi)= I ^dg^Pim^e-^^^'-^'^ 
5 (hi - h{hi-\; gi tjl fj . 



(78) 

Finally, calling respectively A^ , A^ the largest eigen- 
values of the matrices (T n ) p and we obtain that for 
large r 

(79) 

with 



(n n r )£ ~ exp(-r/£p) 
1 



ln(|Al p) /A d | 



(80) 



Notice that all this discussion of stability of the 1- 
rsb solution is not just academic. Indeed simulations per- 
formed on the BM with the distribution P of the fields gij 
generated by population dynamics show that, at a fixed 
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chemical potential /j, > /j,d, the correlation length ^(m) 
increases when to is decreased from m = 1 to 0. In ad- 
dition, the critical length 1/ In 2 is reached at some finite 
value of to, to c < to s , indicating that the description of 
metastable states corresponding to to < m c requires to 
break the replica symmetry beyond one step. 

The limits of the 1-rsb approach will be discussed in 
much more details in the following section devoted to the 
/i = oo limit where it is shown that the [((d — l)oo) x ((d — 
l)oo)] transfer "matrices" (the oo stands for the contin- 
uum range of the fields hi so the matrices are actually 
operators) reduce to finite [(d — l)d x (d — l)d] matrices 
whose eigenvalues can be computed without resorting to 
the population dynamics. 



5 Close-packing limit 

The zero temperature limit of the cavity method, which 
corresponds in lattice glasses to the [i = oo limit, has re- 
ceived particular attention [24], both because of the sim- 
plifications it allows and because of its applications to op- 
timization problems. For lattice glasses, the corresponding 
optimization problem, called the close-packing problem, 
consists in finding, for given lattice and packing constraint, 
the largest achievable particle density. We will obtain the 
solution of this problem as a result of our study. 



5.1 One-step rsb Ansatz 

We first consider the close-packing problem as a limiting 
case of the previous considerations; thus for fi — > oo, the 
single rs equations (3.1) and (12) simplify to 



a = max I 0, 1 — aj + max bj ] , (81) 
j=1 



bo = max I 0, 1 — ^ < 



(82) 



The advantage of this limit is that we can resort to the 
exact Ansatz 



P(a, b) =p e 5{a)5{b) + Pu S(a - l)S(b - 1) 
+ (l- Pe - Pu )5(a-l)S{b). 



(83) 



Note that the simple form obtained here results from 
our appropriate choice of the local fields. Other choices 
may lead to a similar Ansatz but with more than three 
spikes. This minimal number of three is related to the 
three "degrees of freedom" of our model, as appeared 
clearly when we needed three conditional partition func- 
tions. (In general the number of spikes will be the mini- 
mum necessary number of local fields plus one). 

This Ansatz is certainly the only one with integer fields, 
and at this stage it is not obvious why we should not con- 
sider other solutions of Eq. (81-82) with non integer fields. 



However the reason to take integer fields appears clearly 
when working directly at /i = oo. Indeed, the local fields 
have then a simple interpretation in term of the number 
of particles and must therefore be integers. To see why, 
go back to the recursion on rooted trees and note N^ e \ 



N^ u> and the numbers of particles of a (finite) rooted 
tree when its root node i is empty (e), occupied but the 
constraint unsaturated (u) and finally occupied and the 
constraint saturated (s), i.e., the root site has I neighbor- 
ing particles. Considering as before the I = 1 case, we 
have 



N, 



N (e) («) (s) 



k 



1 + 



(e) 



(84) 



= 1 + max V JV< 

U l<]<k 3 P 



e) 



P^jn 



where = raaxi<j<k Nj U \ Obviously, this is nothing 

but the corresponding \i — > oo limit of the Eq. (6-8) with 
SW ~ exp(/zA^ a )), a = e,u,s. The corresponding local 
fields are 



max(7Vf ) ,ivf ) ,ivf ) ) 



7V (e) i 



max 



N 



— N 



r(e) 



(85) 



It is now clear that we can only have a,i,bi € {0, 1} with in 
addition bi < <ij. Moreover, one has a simple interpretation 
of the three spikes. 

Plugging the Ansatz in the general 1-rsb cavity equa- 



tion Eq. (43) and taking y 
parameter we get 



lim^oo fim as breaking 



Pe = Z- 1 (l-p k e -kp k ~ 1 Pu ), 

Pu = z- x p k e ey, 

Z=l + (e«-l)(p k e +kp k e - 1 p u ). 



(86) 
(87) 
(88) 



These equations arc in fact very simple and can be found 
following the principle that for p, — oo a particle must be 
present whenever it is allowed. So in terms of the state of 
the root, the merging of rooted graphs gives: 



e + • • • + e 



u 



(89) 



fe-i 

all other combinations 



which simply means that k empty sites lead to an occupied 
unsatured site, k — 1 empty sites with one occupied un- 
satured leads to an occupied saturated sites and all other 
cases yield an empty site. Now call p e , p u and p s the prob- 
abilities to be respectively in states e, u and s; the three 
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rules translate into three self-consistent equations 

p u oc p k e e v 

Ps oc kp u v k - X e v (90) 
Pe oc 1 - p k e - k Vu p k - x 

which are exactly the 1-rsb equations (86-87) for \i = oo 
when the normalization p e + p u + p s = 1 is taken into 
account. The reweighting factor e v is introduced each time 
the recursion adds a particle since in this case we have a 
"density shift" AN = 1. 

As before, the equilibrium value of the grand potential 
is given by the maximum at y = y s of <f>(y) with here 

-#(y) =ln [1 + (e» - l)( P k e +1 + (k + l)p k ePu )] 

->l±lln[l + (e-y-mi-p e ) 2 -pl)]. 

(91) 

The configurational entropy S(p) as a function of the 
density p is defined by the parametrized curve 

S{y) = y 2 d y cf>( y ), (92) 
p(y) = -d y [y<j>(y)}. (93) 

Note that since d y <j){y s ) = 0, we have for the solution of 
the close-packing problem 

Poo = p(y s ) = -<P(ys)- (94) 

The complexity £ has a maximum at pd = p{y<i) where yd 
is given by d 2 [y(j>(y)] = 0. The curve displays a non con- 
cave part for y < yd which may not have any physical in- 
terpretation; anyway, we will see that this unexpected part 
belongs to a region of the parameter y where the results of 
a 1-rsb calculation are unreliable. In Fig. 6, we present the 
complexity curve for the BM. The close-packing densities 
of various models are presented in Fig. 7. 



5.2 Stability of the 1-rsb Ansatz 

Having derived the <f>(y) function analytically, it is inter- 
esting to compare it with the output of the population 
dynamics algorithm. As a first consistency check, it is ob- 
served that when the population is started on the integer 
spikes, both approaches lead exactly to the same result. 
However, considering the stability of the integer Ansatz 
under population dynamics provides additional features. 

For large y, y > y^ , it is found that even when starting 
with arbitrary fields, the dynamics converges to the ex- 

(2) 

pected distribution on the integers; however for y < y c , 
this Ansatz is found to be unstable and the population 
dynamics converges to a new continuous distribution, cor- 
responding to a greater </>(y), as displayed in Fig. 8. This 
behavior looks puzzling at first sight since on the one hand 
we know that the fields must be integer, and on the other 
hand the cavity method is known to lead to a lower bound 









Pd Pi: 21 Ps\ 



0.57 0.575 0.58 0.585 0.59 

P 

Fig. 6. Complexity curve S(p) obtained from the 1-rsb Ansatz 
for the BM in the close-packing limit p — oo; its slope corre- 
sponds to —y. Four parts of the curve can be distinguished. 
First, a negative part (y > y s ) due to the contribution of un- 

(2) 

typical graphs with few frustrating loops. For p c < p < p a , 
i.e., y^ < y < y s (bold part), the 1-rsb Ansatz is stable and 
p a where S — gives the close-packing density. The complex- 
ity curve for p < p s corresponds to metastable states; it is 
no longer correctly described by the 1-rsb Ansatz for p < p^ 
(y < j/c 2 ^ ) • Finally, we obtain for the fourth part an unphysical 
non concave branch. 



1 


k 


Ps 


Pd 


o {2> 

PC 


1 


2 


0.575742 


0.5703 


0.5739 


1 


3 


0.517288 


0.5097 


0.5159 


1 


4 


0.473384 


0.4646 


0.4728 


1 


5 


0.438382 


0.4288 


frsb 


2 


2 


0.735050 


0.7302 


0.7337 


2 


3 


0.636187 


0.6256 


0.6223 


2 


4 


0.573723 


0.5606 


0.5701 


2 


5 


0.527301 


0.5129 


0.5247 


3 


3 


0.776695 


0.7748 


0.7682 


3 


4 


0.680316 


0.6660 


0.6755 


3 


5 


0.617160 


0.6001 


0.6123 


4 


4 


0.805338 


0.7945 


0.8033 


4 


5 


0.713982 


0.6972 


0.7088 


5 


5 


0.826487 


0.8140 


0.8245 



Fig. 7. Close-packing densities poo = ps for random regular 
graphs of connectivity k + 1 (k = 5 approximates the three di- 
mensional cubic space) where each particle can have no more 
than I neighboring particles. We indicate the dynamical den- 
sity pd where the complexity is maximum; any local algorithm 
trying to determinate p a will stay in the region where p < p d . 
We emphasize however that this value is only a 1-rsb approxi- 
mation (possibly an upper bound) which we have shown to be 
wrong due to the instability of second kind toward further rsb; 
Pc 2 " 1 gives the value of the density where this instability occurs 
and thus provides a lower bound for the correct p d . Note that 
for £ = 1, k = 5 even the equilibrium density p a is not correctly 
described by an 1-rsb Ansatz. 
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Fig. 8. (j>(y) for the BM (£ = 1, k = 2) in the close-packing 
limit fi — oo. The bold line is the result of the 1-rsb Ansatz. 
Its maximum at y a ~ 5.56 gives the close-packing density 
p s = —cj)(y a ) ~ 0.5757. The points with error bars were ob- 
tained with the population dynamics algorithm after 1000 it- 
erations of a population of 10000 fields. We clearly obtain two 
different results for low y. Note that when y increase, error 
bars increases, due to larger and larger reweighting factors, 
making 4>{y) an average dominated by a few large terms only; 
this is why we do not display population dynamics results for 
large y where anyway we know that the 1-rsb Ansatz should 
be recovered. A numerical check of this point is provided by 
studying directly the distribution of the fields. In the inset, we 
show how the fraction of the population within 8 = 0.01 of one 
of the three spikes predicted by the 1-rsb Ansatz evolves with 

(2) 

y. We thus verify that when y > yl ' ~ 5.06, the population 
is entirely on the peaks, even though it was started with an 
arbitrary distribution. 



of the grand potential [25], so given two different solu- 
tions for 4>(y), we must choose the larger. The explanation 
for this contradiction must be that the approximation we 
used, namely the 1-rsb formalism, is not valid. Therefore, 
we expect that further replica symmetry breaking occurs 

(2) 

for the metastable states with y < y c , a situation that 
has been argued to be a generic feature of discontinuous 
spin glasses [26]. 

(2) 

We now present how the exact value y c can be com- 
puted. As when dealing with the liquid instability, two ap- 
proaches are possible. We can either resort to the stability 
of the local fields distribution under the cavity recursion, 
which requires to place oneself in a two-step cavity for- 
malism, or we can stay at the one-step level and consider 
diverging response functions. The fj, — > oo limit of the for- 
malism described in Sect. 4.9 is a bit tricky because cor- 
relations are trivial in the (i — oo limit: (non r ) c = 0, due 
to a total freezing within each state. The stability analy- 
sis based on response functions is still possible, but must 
not rely on susceptibilities; this procedure is presented in 
Appendix B. Here, we will adopt an equivalent approach 
based on the 2rsb formalism, following [26]. 



5.2.1 Two-step replica symmetry breaking cavity method 

To emphasize the generality of the discussion, let us con- 
sider a generic 1-rsb solution at infinite fj, given by a field 
distribution 



P*(h) = J2PaKh - K) = J2PaSa(h) 



(95) 



peaked on d cavity fields, each having (d — 1) components 
that are integer valued [d = 3 for all (k,£) lattice glasses]. 
The p a satisfy a relation analogous to Eq. (86) and (87), 



Pa 



1 

Z 



p bl ...pb k exp[-yAE bl ,... tbk }. (96) 



(&!,...,&*) — >a 



Note that the fields are such that there is no degeneracy, 
i.e., k parents in configurations (&i, . . . , bk) lead to a child 
whose configuration can only be a. 

At the two-step level, not only the configurations are 
grouped into different states, requiring to consider a dis- 
tribution P(h) over these states, but the states are them- 
selves organized into larger clusters, that is groups of states 
sharing some common properties. We therefore need to 
consider a probability distribution Q[P] over the proba- 
bility distributions P(h). It has the following meaning: on 
a given site, the distribution P c (h) of the cavity field inside 
a cluster c must be taken from the distribution Q[P]. 

The 2-rsb cavity equations are obtained by general- 
izing the 1-rsb distribution Af N (w) <~ e Ny(w-u> ) f tne 
number of states with fixed lu to 



J dL}ldLd28(LU—L}l—LU 2 )e 



A/jv(w) 

(97) 

Here ui is decomposed into u>i + u>2 with u>i the grand 
potential of a cluster with respect to a reference u>o, and 
LU2 the grand potential of a state inside the cluster with 
respect to that of the cluster u>\. The hierarchical rsb 
scheme is here reflected by the similarity between the dis- 
tributions of u>i and u>2, (u>i\u>o) ~ exp[A r yi(wi — u>o)] 

and Aftf\(jJ2\vi) ~ cxp[A r y 2 (^2 — Starting from the 
distribution given by Eq. (97) and following the lines of 
the derivation of the 1-rsb cavity equation described in 
Sect. 4.4, we obtain the 2-rsb cavity equation 

^t P ] = | / T\ VP MP3]z[{Pj}] m ' m S[P-P[{P3}]\ 

(98) 

where 



k 

H{Pj}] - j^y } J ft dhjPjihM. h({hj})) 



(99) 



and 



exp(-y 2 Z\£({^})) 

*[{Pj}] = J n^(^ ex P(-^({M))- (ioo) 
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The Parisi parameters yi and y 2 have to be taken such 
that j/i < y 2 . Note that in particular, for yi = and 
V2 = y {y denotes the 1-rsb parameter), the formalism 
describes a non-factorized 1-rsb solution. 

It is essential to understand how the 1-rsb Ansatz must 
be written in this 2-rsb formalism. Two scenarios are in- 
deed possible. Either the 1-rsb states coincide with the 
2-rsb states and there is just one trivial 2-rsb cluster, or 
the 1-rsb states coincide with the 2-rsb clusters and 2-rsb 
states reduce to single configurations. Within the first sce- 
nario, the one-step corresponds to Q = S[P — P*} where 
P* = ^2 a Pa.da is the one-step probability distribution of 
Eq. (95) while within the second scenario, the one-step 
corresponds to Q[P] = J2 a P^[ p - 

Depending on which case we consider, we can have two 
possible kinds of instabilities, as first noted by Montanari 
and Ricci-Tersenghi [26]. As in Sect. 4.9, in the first case, 
the states gather into different clusters, while in the second 
case new states appear within the old states which there- 
fore become clusters. A pictorial view is given in Fig. 3. 

5.2.2 Instability of the first kind: aggregation of states 

The instability of first kind can be studied by considering 
an Ansatz of the form 

Q[P] = f[P - P*] (101) 

where / is a functional with support around the null func- 
tion. The instability is given by the eigenvalue of largest 
modulus A\(y) of the Jacobian matrix associated with 
Eq. (96). Here again, as we deal with random graphs we 
ignore the modulation instability k\Ai\ > 1 and focus on 
the glass instability \/fc|yli| > 1. Different cases are ob- 
served as we vary the parameters I and k in our lattice 
glass models. In some cases the instability is absent and 
appears only asymptotically, i.e., we have v^yli^)! < 1 
for all y but \/fc|^li(y)| — > 1 as y — > oo; this happens on 
low connectivity graphs, e.g. for k — 2,3 when I = 1. 
At higher connectivities, we can define a critical y^ such 
that VklA^l > 1 for y > yi 1} . Then we have to de- 
termine the relative position of y^ with respect to y s 
giving the maximum of </>(?/)• Fig- 9 shows how y s and 
y^ evolve with the connectivity k + 1 for the case 1=1. 
When y s < as it is found for k = 4, 5 for £ = 1, the 
positive part of the complexity curve is unaffected and we 
can rely on our 1-rsb description for typical graphs. How- 
ever, if y s > y^c \ as we find when 6 < k < 25, the 1-rsb 
treatment is not stable, and one should develop a higher 
order rsb formalism. 



5.2.3 Instability of the second kind: fragmentation of states 

To study the instability of second kind, we consider an 
Ansatz of the form 

Q[P} = J2Pafa[P-S a ] (102) 

a 



where the f a have support around the null function. Since 
P is necessarily a combination of the S e , the argument of 
f a can be written as SP = P — 8 a = J2 e ^a e e(^e — 5 a ) and 
we evaluate the widening of the f a by computing (5P) a = 
J T>SPf a (5P) to obtain the following relations: 

P a IU u \ 

( 103) 

E e-v* AE ^-^(e c ) bl . 

c^6i,(c,62,...,6fc)— >e 

Since we consider a local instability of the 1-rsb solution 
we take yi = y 2 = y- Noting X a ^ e = (e e ) a it reads 

PaK^e £ Phi ■ ■ ■ Pb k 

(6i,...,6ji)— »o 

£ e-^- >>\ b . , n . 

j,c^bj,(bi, ...,c,...,6 fe )^e 

(104) 

This notation emphasizes the relation with the point of 
view based on response functions developed in Appendix 
B: p(e\a) = p a ^a^e is the probability to replace a with e 
by changing only one parent. Indeed, for /j, — > oo, adding a 
small pcrturbating field is equivalent to changing a config- 
uration a. The equation tells us how such a change prop- 
agates from a site to its neighbor. 

As expected from Sect. 4.9, the instability is described 
by a [d(d — 1)] x [d(d — 1)] transfer matrix T ac i,bc, ( £ e)a = 
k*llb^cTae,bc( € c)b- If ^2 is the eigenvalue of T of largest 
modulus, the one-step solution is stable provided kA 2 < 1 
(A2 is positive, due to Perron- Frobenius theorem) ; here we 
need not square the eigenvalues since the first moments 
(e e ) does not vanish due to the positivity of the A a ^ e 
(since A 2 is positive the instability detected by kA 2 = 1 
is compatible with the underlying random lattice). Note 
that the same principles are straightforwardly extendable 
to cases with fluctuating connectivity and/or quenched 
disorder; we would have to average over the different pos- 
sible transfer matrices. 

We present in Appendix C a simple method to ob- 
tain the transfer matrix T and focus here on the results. 
Generically, it is found that kA 2 < 1 for large enough y, 

(2) (2) 

y > y c ■ For instance for the BM, we obtain y y c ' ~ 5.061. 
This instability is precisely the one detected by the pop- 
ulation dynamics algorithm. To underline this point, we 
determine the fraction of the population within 6 = .01 of 
one the d = 3 expected delta peaks after 1000 iterations. 
At y = 5.06, this fraction is found to represent 99.2 % of 
the total population composed of 50000 individuals, while 
for y = 5.07 (and all larger values) it is found to be 100 % 
as shown in the inset of Fig. 8. We can thus conclude that 
the onset of non integer fields in the population dynamics 
is a clear sign that the replica symmetry must be broken 
beyond one step. 

(2) 

When comparing y c with y s , the value of y for which 
(f>(y) taken with the 1-rsb Ansatz is maximum, we find 
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Fig. 9. For I = 1 and various k, values of the parameters y s , 
yiP and yi 2 - 1 giving respectively the equilibrium thermodynam- 
ics, the instability of the first kind (for y > yc) and of the 
second kind (for y < y^). Only if y^ < y s < y^ is the close- 
packing limit correctly described by the 1-rsb Ansatz. That is 
the case for k < 4 and k > 26. The lines are only guides for 
the eyes (we omitted y c ^ ^ 7.0 for k = 4). 



(2) 

either y c < Us indicating that the close-packing is 1-rsb 

(2) 

or > y s suggesting that it is full rsb. For a given t, 
we obtain that y^ < y s at low connectivity (k < k\, with 
fci = 4,8, 10 for respectively £ — 1,2,3 ). But things are 
not exactly that simple since for instance when I = 1 and 

(2) 



k > ki = 19, y c < y s again 



5.3 Nature of the close-packing 

We now summarize the possible nature of the close-packing 
by taking into account both kinds of instabilities. Sev- 
eral situations can occur, corresponding to the different 
relative positions of y s , the maximum of <p{y) giving the 
equilibrium properties with metastable states described 
by y < y s , y^ indicating an instability of the first kind 
for y > y c X \ and y c 2 ^ indicating an instability of the sec- 

(2) 

ond kind for y < y c (note the different directions). All 
possible combinations seem to be already contained in the 
£ = 1 model so we will detail it for connectivity ranging 
from k = 2 to k = 26; illustrations are provided by Fig. 9 
and 10. 

Low connectivities k — 2,3,4 make the best glass mod- 
els since yjp < y s < y^, which means that the equilib- 
rium state is 1-rsb and only less dense metastable states 
are frsb [case (a) of Fig. 10], a behavior similar to p-spin 

models [26]. When k = 5, y s < y c < yi so the equilib- 
rium state is frsb; however there exist some (but not all) 
untypical graphs with a 1-rsb phase [case (b) of Fig. 10] . 
When 6 < k < 24 the 1-rsb Ansatz is never stable since 
J/c 1 ^ < y^] however the situation improves somewhat as 
k increases, with y 1 ^ < y s < y^ for 6 < fc < 18 ("dou- 




(b) 



(c) 



(d) 




Fig. 10. Part of the complexity curve (bold line) correctly 
described by the 1-rsb Ansatz. For I = 1, the case (a) is realized 
for k = 2, 3, 4; the models with k = 2 (BM) and k = 3 have 
in fact no instability of the first kind, i.e., the whole p > 
range is stable. Case (b), realized for k — 5, corresponds to a 
frsb close-packing, as well as case (c) and (d). In case (c) no 
part of the l-rsb complexity curve is correct; this happens for 
6 < k < 18, the same conclusion holding for 19 < k < 24 where 
Pc 1 ' < p s < p^? [case (c'), not represented here]. In case (d), 
realized for k = 25, some of the metastable states are indeed 
described by the 1-rsb Ansatz even if the close-packing is frsb. 



for 19 < k < 24 [case (c) of Fig. 10]. The case k = 25 
is particularly interesting and justifies our study up to 

this large connectivity; indeed, we find y^ < y^ < y s , 
i.e., the equilibrium state is frsb as well as the densest 
metastable states but there also exist metastable states 
with intermediate density which are well described by the 
1-rsb approach [case (d) of Fig. 10]. Such a strange fea- 
ture disappears for k = 26 (and presumably for all higher 
k) where we recover the low connectivity situation with 



bly unstable") and y^ < yf < y s ("simply unstable") 



(2) 



yf ] <y s < yc 1] - 



6 Dynamical and kinetic transitions 

As explained in Sect. 4.8, our lattice glass models undergo 
a dynamical transition at /id where the phase space breaks 
into exponentially many metastable states. Here "dynam- 
ical" means that the transition does not affect equilib- 
rium properties, i.e., the grand potential is analytic at jid ■ 
However, despite its designation, the dynamical transition 
refers to a static property, in the sense that it describes 
a morphological change of the free energy landscape and 
makes no reference to a particular kinetic rule. Kinetic 
transitions (also called "dynamical arrests" [27]) are as- 
sociated with the freezing of some degrees of freedom of 
the system under specific local kinetic rules, and form the 
central concept of kinetically constrained models [2] . 

We now show that choosing the most natural local 
dynamical rules leads to several properties: (1) A kinetic 
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transition at which a finite fraction of particles is blocked 
takes place beyond a certain chemical potential [ik- (2) 
The kinetic freezing at [ik is a distinct phenomenon from 
the ergodicity breaking occurring at [id- In particular it 
may happen that [ik < fJ-d so that the dynamical transition 
at Hd is prevented by a kinetic arrest at [ik- 

The easiest and more natural dynamical rule for our 
models (the one that, by the way, one can use naturally 
to simulate it) is the standard Monte Carlo, i.e., at each 
time step we pick at random a particle on a site % and 
one of its neighboring site, j, and we move the particle 
onto j provided that the geometrical constraint remains 
satisfied after the move (site j has no more than I occupied 
neighboring site in addition to site i). Furthermore, for 
simplicity, we focus on the simple I = 1 model and instead 
of computing the exact kinetic transition [ik as done in 
Ref. [4] for the Kob- Andersen model on the Bethe lattice, 
we present a simpler calculation consisting in finding the 
onset of a particular blocked structure. This provides an 
upper bound fib > [ik for the kinetic transition, so that 
obtaining [id > Hb will be enough to prove that [id > [ik 
can occur. 

The blocked structure we consider is a percolating struc- 
ture made of two kinds of blocked particles, either unsatu- 
rated or saturated (see Fig. 11). Rooted trees with an un- 
saturated and saturated blocked particle on the root are 
associated with the partition functions S( bu ) and E^ bs ^ re- 
spectively. The relation between the partition functions at 
different generations g reads 

-3+2 - e 



fc-1 



"3+3 - Ke 



fc-1 



(.9 -> g + 1) 



fc-i 



(105) 



where E^> = E^ + E^ , and E = E^ + E^ + E^ is 
the total partition function (with notations of Sect. 3.1). 
As illustrated in Fig. 11, the equation for E g b ^l includes all 
configurations at generation g [term (E g ) k ] except those 

where no particle is blocked \(yE g — Eg b ^ ] and those 

where only one particle is blocked [kE^ [e 9 — Eg b ^ ]; 

in this case, the presence of two blocked particles on level g 
is needed to exclude the possibility that a blocked particle 
on level g could go at level g + 1 and then to a differ- 
ent site at level g. The equation for E^l has a similar 
interpretation. 

In the liquid phase where the partition functions be- 
come independent of the generation <?, we can consider 
the probability for a particle on the root to belong to the 




# bs 



Fig. 11. Illustration of the blocked structure under consider- 
ation, with empty sites represented in white, sites occupied by 
a blocked particle in black and unspecified sites in grey. The 
structure is made of two kinds of blocked particles. The left 
figure shows how unsaturated blocked particles (bu) at gen- 
eration g + 2 are prevented from moving to generation g + 1 
provided sites at generation g + 1 have at least two blocked 
neighboring particles b (b meaning either bu or bs). The right 
figure shows how saturated blocked particles (6s) at generation 
g + 3 are similarly blocked by blocked particles at generation 
g + 1 and g. In this case a branch can be blocked either by a 
single 6s particle, or by two blocked particles. 



blocked structure, 



Pb = — 



rJ(«) + E( s ) ' 



(106) 



From Eq. (105) and (6-8), it is found to satisfy the self- 
consistent equation 



Pb = V 
+9 



1 - (1 - (Pbf - k(r[ (l - (1 - CPbf (107) 

-fcCP6(l-CP5) fc_1 ) fc (l-CP6) fc_1 



2fc-l 



where £, rj and 9 depend on fi through the relations ( = 
I _ e-f"".-, r) = ( e ^ii<. - l)/( e ^ a "o - 1) and 9 = (e^ - 
e nbuqy (g^anq _ 1). At a given connectivity k we calculate 
[ib as the lowest [i such that Eq. (107) has a solution pb € 
]0, 1]. For k < 5 we obtain [ib > [id, he., the blocked 
structure considered does not appear in the liquid phase. 
However, for k = 6 we find [ib — 1.0 which is lower than 
[id — 2.7; we can thus conclude for sure that [ik < [id, ke., 
the kinetic freezing occurs while the system is still in its 
liquid phase. 

Let us conclude this section with some remarks on the 
difference between the dynamical transition obtained from 
the cavity (or replica) computation and the kinetic freez- 
ing that we studied in this section. The kinetic freezing 
transition is clearly dependent on the dynamics and is re- 
lated to the fact that after a certain density there are no 
more paths to go from one part of the configuration space 
to another. Of course, the allowed paths depend on the 
local dynamics that has been chosen for the model. If one 
increases the scale on which particles can move, for exam- 
ple allowing particles jumps on next nearest neighbor or 
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further, then the density pk, at which the kinetic freezing 
takes place, is expected to increase. The thermodynamic 
limit, which does not play an important role for kinetic 
freezing transitions, is instead very important for the dy- 
namical transitions due to a change of the morphology 
of the free energy landscape. In the latter case when the 
number of sites becomes very large some bottlenecks in 
the configuration space (that has to be used by any local 
dynamics) shrinks. So, at the end, the configuration space 
breaks up in different ergodic components. 



7 Discussion and conclusions 

Some of the qualitative properties of our models such as 
the presence of a liquid phase at low fi are easily under- 
stood. Since our work concerns random graphs with no 
local disorder, the equilibrium liquid phase is character- 
ized by its homogeneity: the probability that a site is oc- 
cupied is site-independent. Within the cavity framework, 
this is reflected by a unique value of the cavity fields which 
describe in fact probabilities of particle occupation. Due 
to the tree-like structure of the lattices we consider, such 
probabilities can be computed recursively. However such 
an approach is physically justified only if the iterations 
converge to a unique fixed point regardless of the initial 
conditions. 

When fi is increased, the liquid fixed point becomes un- 
stable and so some kind of equilibrium glassy phase must 
be present. A clear evidence of the inadequacy of the liq- 
uid solution (i.e., of the Bethe-Peierls approximation) is 
the negative value of the liquid's entropy. In some cases, 
crystal phases can be constructed and we found these to 
be thermodynamically favored beyond a melting point fi m ; 
however such ordered phases occur only on trees with care- 
fully chosen boundary conditions. By focusing on random 
graphs, we exclude this possibility, leaving as the only al- 
ternative the existence of an equilibrium glass phase. 

The nature of the liquid to glass transition is discon- 
tinuous and of the one-step rsb type. This means that the 
phase space first clusterizes into exponentially many non- 
ergodic components at some fid, before a true equilibrium 
transition occurs at fi s > fid- Both transitions are ther- 
modynamical in nature but the transition at fid is called 
a "dynamical transition" to emphasize the fact that equi- 
librium properties are still given by the liquid solution 
when fid < fi < fi s , and only non-equilibrium properties 
are affected in this range. More insights into these non- 
equilibrium effects where obtained by computing the con- 
figurational entropy S{uj) which gives the number Af(u>) of 
(metastable) clusters of configurations with a given grand 
potential u), through the relation Af(u>) = exp[NS(ui)] 
(with N being the number of sites of the graph). 

Investigating the stability of the 1-rsb solution, we 
found that while the equilibrium configurations may be 
of 1-rsb type i.e., organized into distinct clusters, some 
metastable states can be associated with full replica sym- 
metry breaking (frsb), i.e., configurations organized into 
clusters themselves organized into smaller clusters, and so 



on. Our study of the fi = oo case furthermore demon- 
strated that, depending on the model, the equilibrium 
state itself can exhibit frsb, and many different phase 
space structures are possible. 

Finally, we showed that when using standard Monte 
Carlo dynamical rule a kinetic freezing can occur at fik < 
fid, thus preventing the dynamical transition which is due 
to a change of the free energy landscape. In a general 
model, the freezing will therefore occur at fif — min(pk , fid) 
and will be of a different nature, purely kinetic or truly 
thermodynamic, depending on which alternative is real- 
ized. The overall picture for lattice glass models given here 
is quite rich. A natural question is how small loops modify 
(if at all) this picture; we expect only small quantitative 
differences. A more fundamental and intriguing issue is 
therefore what features survive when going to finite di- 
mensional lattices where the large scale Euclidean struc- 
ture could potentially change even the qualitative aspects. 

The last point we would like to emphasize here is that 
even if our conclusions do not necessarily apply to real 
vitreous materials, our results should have an impact in 
optimization problems where random graphs arise natu- 
rally. Indeed, as a by-product of our study, we have deter- 
mined the maximum densities and the phase space struc- 
ture of different close-packing problems. Our family of 
models provides a generalization to I > 1 of the hard- 
core model t = (called the vertex cover problem in com- 
puter science and the largest independent set problem in 
the mathematics literature) with, in some cases, 1-rsb fea- 
tures much simpler than the frsb structure of the simplest 
(£ = 0) hard-core model. A full understanding of these 
models could therefore constitute a first step toward a 
mathematical resolution of this fundamental but still un- 
solved problem [28] . 
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Appendix A 



In this appendix, we present the generalization of the for- 
mulae given for I = 1 in the core of the paper to an ar- 
bitrary £ < k, where t is the maximum number of neigh- 
boring particles a given particle can have. It is straight- 
forward to realize that here again, only three degrees of 
freedom are relevant (d = 3) so that we also need only 
two-component local fields and the situation is very simi- 
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lar to I = 1 . The recursion relations are now 
a =a(ai,bi, . . . ,a k ,b k ) 

=Iln [i + e Mi-Ej =1 a,) 
u L 

!+E E fl(^-i) 

p=l l<ii<---<i p <fc r=l 

6 =6(01,61, • • • ,a k ,b k ) 

=-ln h+e^-S^a,) 
u L 



(108) 



(109) 



i-i 



!+e e n^- 1 ) 

p=l l<ii<---<i p <fc r=l 

In the jjL — > oo limit this becomes 



= max I 0, 1 — > a,- + 

/ 

6 = max 0,1— > a, + 

\ ^-f J i<ii 



max &j 

r=l 
(-1 

max > b; 

r=l 



(110) 



The Ansatz on integer peaks has exactly the same struc- 
ture as for i = 1, 



P(a, 6) =p e 5(a)S(b) + Pu S(a - l)5(b - 1) 
+ (l-Pe-Pu)6(a-l)8(b) 



(111) 



but of course the relations between the p a need to be ex- 
tended: 

S=l + (e»-l)£(*V 

9=0 



(112) 



Next we can write the corresponding (j)(y) as 



-#(y)=ln(l + ( e «-l)W* 



fc + 1 



ln(l + (e-"-l)[(l-po) 2 -p?]). 



(113) 



The study of the stability follows the same principle as 
for I = 1, the only difficulty being to correctly collect all 
the combinatorial factors. The matrix for the instability 



of the second kind is 



T = 















^01,20 


^01, 21 








*02, 10 





t<)2, 20 


t02,21 





tl0,02 





tl0,12 








£l2,01 



















(k- l)e- y * 



















I 












with 



(114) 



P k - q - 2 p q u p s e\ 



9=0 \ y / 

lE^zJV^rw* 

9=0 \ q ' 



toi, 20 

^01,21 
*02,10 



1 /fc — l\ 

*02,20 = -^( ; J^-'- 2 p^ s e«, 

<te, a i = |*(Jl 1 1 )rf- , -V«-W y , 



(115) 



t 



10,02 



10,12 



1-12,01 



Z ^o 



g=0 

<j=0 vy/ 



where p s = l-p e -p u and Z = l+(e^-l) J2 q=0 QP q J t e~ q '■ 

Appendix B 

In this appendix we shall show how the stability criterion 
at fi = oo derived in the text can be obtained within a 
response function formalism. 

First let us note that as far as the fragmentation of 
the cluster is concerned we can merely take the = oo 
of the formalism introduced in Sect. 4.9. Indeed in the 
close-packing limit there is a total freezing within each 
state. Hence, it is relatively easy to obtain the correlation 
length £2- The computation of (nitij) can be done using 
the transfer matrix: 

T A ,B = Pa 2-- - Pa " CX Vhy AE (A,a 2 ,...,a k )}, 

(A,a 2 ...,a k )^B 

(116) 



(n n r ) = 



pt T rp> 



(117) 



where Pq and Pq are suitable vectors to impose that there 
is a particle in and r. The other vectors P* and P* read 
with the bracket notation: 

(P*|a)=p a , (a\Pi)= ]T Pb . (118) 

b compatible with a 
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What is important to notice is that they are respec- 
tively the left and right eigenvectors of T with largest 
eigenvalues, Z, in module. The properties of being eigen- 
vectors can be checked but it is fundamentally due to the 
fact that p a is the self-consistent solution of the itera- 
tive cavity equations which can be written going from left 
to right or from right to left on the ID chain. The fact 
that they correspond to the largest eigenvalue in module 
is just due to the Perron-Frobenius theorem (T is irre- 
ducible with positive coefficient and the two eigenvectors 
have all the components with the same sign). 

Gcnerically T can be written in terms of its right and 
left eigenvectors: 

T = ^A a \v a >< w a \, < wp\v a >= 6 a ,p, (119) 

a 

where the a = 1 corresponds to P* , P* (note however that 
P* and PI are not normalized between one another) and 
At = Z. 

Plugging this expression into the equations for the cor- 
relation function one obtains generically that the correla- 
tion function reads: 

(W) =C + Y,Di (^) (120) 

i=2 ^ 1 ' 

where C and Di are constants that can be expressed in 
terms of the scalar product of the eigenvectors. In par- 
ticular it can be shown that C — p 2 as it should. Thus, 
the stability with respect to fragmentation of the cluster 

is governed by the condition k (^jf'j < 1- 

Finally, we want to show that this condition is the same 
as the one given in Sect. 5.2.2. Let us define the matrix T: 

Ta,b = \ [t a , b -PbJ2 T A,cj ( 121 ) 

This matrix has eigenvalues:0, . . . , One can check 

this in different ways. For example TrT' = Y^i=2 (yfr) ^ or 
any I, so the eigenvalues are indeed the ones predicted. In 
terms of T the stability criterion is that k times the square 
of the largest eigenvalue of T has to be less than one. This 
is exactly the criterion obtained in the text by studying 
the stability of the iterative Eq. (86) for the 1-rsb glass 
(where the normalization Z is considered as a function of 
the p a ). 

The second instability, the fragmentation of states, 
cannot be obtained taking simply the u — > oo limit of 
the formalism described in Sect. 4.9. Because of the com- 
plete freezing within a state in the close packing limit the 
connected correlation functions are trivially zero. This, of 
course, does not mean that there is no way to identify a 
correlation length. Actually, one has just to resort to a 
definition of correlation length in terms of response func- 
tions. Let us remark that the correlation function (non r ) c 
can be written as: 

PoPo (Mo - (n r ) e ) (122) 



where Pq and p% are respectively the probability that the 
site is occupied or empty and the averages (-) and (-) e 
are the averages conditioned on the events n = 1 or 
no = 0. At finite fi computing (non r ) c or the response 
function ((n r } — (n r ) e ) contains the same amount of in- 
formation as far as the correlation length is concerned. 
However at infinite /i there is total freezing within a state. 
As a consequence the product PqPq equals zero and one is 
forced to extract the correlation length from the response 
functions. Thus one has to study the change in probabil- 
ity for the occupation variable at site r when the value of 
occupation variable in is changed. This gives rise natu- 
rally to the eigenvalue problem discussed in Section 5.2.3. 
In that case the matrix entering in the eigenvalue prob- 
lem tell us how the effect of replacing a with e by changing 
only one parent propagates from a site to its neighbor. 



Appendix C 

We present here on the specific case of the lattice glass 
with 1 = 1 how Eq. (104) allows for a straightforward 
study of the 1-rsb stability. We start with the composition 
rules Eq. (89). With the notations of this section, we have 
d = 3 and the configurations are a € {e,w, s}. Let us 
see for instance how the matrix elements corresponding to 
\ s ^ e are computed. We have to ask how the configuration 
of one parent on the left hand-side of 

u + e-\ \-e—>s (123) 

must be changed if we want to obtain a configuration e 
instead of s. Begin with the parent u: according to the 
composition rules, if we replace it by e it leads to a child 
u so we do not retain this possibility; instead, if we replace 
it by s we indeed obtain e. If we now keep u but try to 
"flip" one of the e's, we see that both replacements by u or 
s lead to a child e as we wish. So we have a total of three 
cases to take into account, a situation that we summarize 
by writing the following "reaction" rules: 

(e|s) <— (s\u) + ej he , 

fe-i 

(e\s)<-u+(u\e) + e + -y + e , 
k-2 

(els) < — u -\- (s\e) + e + ■ ■ ■ + e . 

* v ' 

k-2 

Now recall that the initial "reaction" formula (123) rep- 
resents in fact the equation 

Ps = Z^kpJT 1 ^ (125) 

where Z is a normalization constant, kpup^," 1 is the prob- 
ability to choose parents (u, e, . . . , e) when each configura- 
tion a is taken with its probability p a , and the rewcighting 
term e v is added each time the generated child is u or s, 
i.e., corresponds to a new particle. Generalizing slightly 
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the same principle, just by looking at the reaction equa- Now consider n replicated systems for the same ran- 
tions (124), we write dom graph and define the function ^>(t): 



fc-2 



p(e\s) =Z 1 kp(s\u)p k e 1 + Z 1 k(k - l)p u p(u\e)p 

+ Z-'k(k l) Pu p{ s \e)p^. *(T) = ^< II eMh^ + g^r*) > (132) 

(126) a=1 

Here we put no reweighting since we generate an empty where the overline means the average over the disorder 
site e. Using the equations for the p a , it can be simplified ( or ' analogously, one can average over sites for a fixed 
somewhat, and when expressed with the \ a ^ e defined by random graph), N is a normalization constant such that 
p{e\a) =p a X a ^e, it becomes £ T = 1 and r = (n, . . . ,r n ). 

For the same reasons as discussed in the text in the 
\ s ^ e = [\ u ^s + (k — l)\ e ^u + (k — l)\ e ^s]e~ y ■ (127) case of the cavity method, the grand potential per site is 

given by: 

2w = lim \ (-(k + l)G<£i + 2G&>) , (133) 

A_ =*A_, G H = ln E iKrWrWr, t>), (134) 

A -A 

fe+1 

\ e ^ u =Z- 1 kp s p k - 2 \ s ^ e e\ G («) = ln J2 J] 1>(Ti)C({Ti}, n ) (135) 

A e ^ s =z- 1 kp s p k - 2 \ s ^ ue y + z~ x k(k - i)p 2 u p k e - 3 \ u ^ ee y *■«,»„ i=i 

+ Z k(k-l)p u p s p e 3 X s ^ e e v . where C(t,t') and C({Ti},n ) enforce the geometrical 

(128) 

constraint. 

These six linear relations define the matrix T; we then ^Differentiating w with respect to V(r) we get the self- 
compute its eigenvalue A 2 as a function of y. consis ent equation 



The five other equations are obtained by following the 
same lines, 



Appendix D 



V(t) = M J2 IT <P(n)C({ n }, n ) (136) 



In this appendix we would like to show how the results , , , iuj.ri.ir 

. . , . . . . . . where the prime over L means that, for each value of r, 

obtained m the text with the cavity method can be re- we gum Qm oyer thoge ^ nQ ^ iye rige tQ T 

produced with replicas. For an introduction to the replica „ , ,. ,, s . , , 

n r , , . ±1 c ., Once the distribution w(t ) is known all the observables 

method tor random graphs with finite connectivity sec tor , , „ . , , u j i 

example [30] Can computed, tor instance the density equals: 

First, we introduce a function to encode in a compact k+1 

way the conditional probability that the central site i for a p = lim — n a TT iP(t-)C({t-} n ) (137) 

rooted sub-tree is empty (e), occupied but the constraint «->o n ^ °! j ' * 
unsaturated (u) and finally occupied and the constraint 

saturated (s), i.e., the root site has I neighboring particles. xhc mQgt ral u symmetric solution can be 

We use the function /(t) — Ccxp(/iT + gr ) with r — written' 
— 1,0,1 where C,h,g are defined in such a way that 

-( e ) if \ f av.j vth \ TT exp(hT a + gr 2 ) 

/(o) - „ (e) x „ w , ^)=y^^^)n 1+exp(9+fc)+expa ,_ fc) 

■=» + + ~i (138) 

where P is a normalized probability distribution. Because 

fi 1 ) = (e) ( M ) ( a ) ■ ( 129 ) of the homogeneity of the random graph with fixed con- 

"* ncctivity, the rs solution is particularly simple, 

/(-l) = 1 - /(0) - /(l) = ~* (u) . P(h,g) = 8(h-h^)8{g-g IUl ). (139) 

The relation with the fields a, b can be easily established: Using the relations (130-131) between h, g and a, b one 

can rederive all the liquid properties (density, free energy, 

exp(— ua) = - (130) entropy, . . . ) obtained in the main part of the paper with 

1 + cxp(g - h) + exp(g + h) ' the cavity method. 

1 Because of the homogeneity of the random graph, the 

ex Pl A* ) — ^ — cxp(g — h) ' 1-rsb solution, in which replicas are divided in n/m groups 
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(noted c) of m replicas each, also becomes simple and 
reads 



V(r) 



n/m 

n 



dhdgP(h, g) 



grl) 



[1 + cxp(.g + h) + exp(g 



- h)] m 
(140) 

Plugging this expression into the general equation for i/j(t) 
and changing variables from h, g to a, b we get back to the 
equation on P(a,b), Eq. (43). Moreover using this 1-rsb 
expression of iP{t) one can easily obtain all the observables 
like the density, the free energy and the complexity. 
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